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Abstract 

In this paper, we apply the method of approximate transformation groups proposed by Baikov, Gaziziv 
and Ibragimov [I] [2], to compute the first-order approximate symmetry for the Gardner equations with 
the small parameters. We compute the optimal system and analyze some invariant solutions of These 
types of equations. Particularly, general forms of approximately Galilean-invariant solutions have been 
computed. 
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1 Introduction 

The Kortewege-de Vries(KdV) equation is a mathematical model to describe weakly nonlinear long waves. 
There are many different variations of this PDE, But its canonical form is, 

Ut - 6UU^ + U^a:x ^ 0- (1-1) 



The KdV equation (1.1 1 owes its name to the famous paper of Diederik Kortewedg and Hendrik de Vries, 
But much of the significant work on the KdV equation was initiated by the publication of several papers 
of Gardner et al. (1967-1974). The remarkable discovery of Gardner et al. (1967) that KdV equation was 
integrable through an inverse scatering transform marked the begining of soliton theory. 
In 1968 Miura [3 HO] introduce "Miura transformation", 

u = v^+v^, (1.2) 

to determinean an infinite number of conservation laws. If we put v = + which e is an arbitrary real 
parameter, then Miura transformation becomes 

u — -e^^ + w + eWx + e^w^ . 

But, since any arbitrary constant is a trivial solution of KdV equation, it may be removed by a Galilean 
transformation, so we just consider "Gardner transfosmation", means, 

u — w + eWx + e^iiJ^ , 

substituting above transformation in KdV equation, observe that w satisfies in "Gardner equation", 

wt- Q{w + e^w'^)wx + Wxxx (1-3) 

for all e. 

Obviously, if e = 0, Gardner equation becomes the KdV equation. Next, we observe that the for more 
detailed exposition and references refer to [31 [S] . 
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These two nonlinear equations are integrable with inverse scattering method. But in this paper, we 
analyse them with a method which is introduce by Baikov, Gazizov and Ibragimov PQH]. This method 
which is known as "approximate symmetry" is a combination of Lie group theory and perturbations. There 
is a second method which is also known as "approximate symmetry" due to Fushchich and Shtelen [3] and 
later followed by Euler et al OE]. For a comparison of these two methods, we refer the interested reader to 
the papers [12l [13] . 



2 Notations and Definitions 

Before continuing we need to present some definitions and theorems of the book of Ibragimov and Kovalev 
If a function f{x,s) satisfies the condition 

hm ^ = 0, 

e-i-O eP 

it is written f{x,e) = o(eP) and / is said to be of order less than e^. 

If /(x, e) — g{x, e) — o(e^), the functions / and g are said to be approximately equal (with an error o(eP)) 
and written as f{x,s) = g{x,e) + o(£^), or, briefiy f ~ g when there is no ambiguity. The approximate 
equality defines an equivalence relation, and we join functions into equivalence classes by letting f(x, e) and 
g{x, e) to be members of the same class if and only if / 5. 

Given a function f(x,e), let 

fo{x)+efiix) + ... + ePfpix), (2.1) 

be the approximating polynomial of degree p in e obtained via the Taylor series expansion of f{x, e) in 
powers of e about e = 0. Then any function g ~ f (in particular, the function / itself) has the form 

g(x, e) = h{x) + eh{x)... + f^{x) + o(eP). 



Consequently the expression (2.1) is called a canonical representative of the equivalence class of functions 
containing / Thus, the equivalence class of functions g(x,e) ~ f{x,e) is determined by the ordered set of 
p+l functions /o(a;), /i(x), /p(x). In the theory of approximate transformation groups, one considers 
ordered sets of smooth vector-functions depending on x's and a group parameter a: 

fo{x,a),fi{x,a),...,fp{x,a), 

with coordinates 

foi^,a),fl{x,a),...,Pp{x,a), i = l,...,n. 
Let us define the one-parameter family G of approximate transformations 

x^«/^(a;,a) + e/^(x,a) + ... + ef/;(a;,a), z = l,...,n (2.2) 

of points X = (x^, ...,a;") £ M" into points x — (x^, ...,x") G M" as the class of invertible transformations 

x = /(x,a,£), (2.3) 

with vector-functions / = (/^, /") such that 

f (x,a,£) w r^{x,a) + £fl{x,a) -|- ... -f eP/^(x, a), i = l,...,n. 

Here a is a real parameter, and the following condition is imposed: 

/(x,0,e) w X. 



The set of transformations (2.2) is called a one-parameter approximate transformation group if 

/(/(x, a, e), &, e) « /(x, a + h,e) 



for all transformations (2.3). 
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unlike the classical Lie group theory, / does not necessarily denote the same function at each occurrence. 
It can be replaced by any function g ~ f. 

Let G be a one-parameter approximate transformation group: 



« f{z, a, s) = faiz, a) + efl{z, a), i ^ 1, N. 



An approximate equation 



F{z,e) = Fo{z)+sFi{z)^0 
is said to be approximately invariant with respect to G, or admits G if 

F{z,e)^ Fif{z,a,e),e)=ois) 



(2.4) 
(2.5) 



whenever z = {z^ , z^) satisfies Eq.(2.5). 

If z = (a:, u, U(fc)), then (2.5) becomes an approximate differential equation of order k, and G is 

an approximate symmetry group of the differential equation. 



Theorem 2.1. Equation (2.5 1 is approximately invariant under the approximate transformation group (2.4) 
with the generator 

X = Xo + eXi=Coiz)g^+ea{z)^, (2.6) 



dz 



if and only if 



X'-''^F{z,s) 



F«0 



o(e), 



X^„'^Fo{z)+e[Xl'"Foiz)+X'„'"F,iz) 



(fc). 



Die 



(2.7) 



In which k is order of equation , and X'^^^ is fc*'' order prolongation of X. The operator (2.6) satisfying 
Eq. (2.7 1 is called an infinitesimal approximate symmetry of, or an approximate operator admitted by Eq. 
(2.5). Accordingly, Eq. (2.7 1 is termed the determining equation for approximate symmetries. 



Remark 2.2. The determining equation (2.7) can be written as follows: 

xl^^Foiz) = A(z)Fo(z), 



(2.8) 
(2.9) 



The factor \{z) is determined by Eq. (2.8) and then substituted in Eq. (2.9). The latter equation must hold 
for all solutions of Fo(z) = 0. 

Comparing Eq. ( |2.8[ ) with the determining equation of exact symmetries, we obtain the following state- 
ment. 



Theorem 2.3. If Eq. (|2.5|) admits an approximate tramformation group with the generator X = Xq + eXi, 

d 



where Xq ^ 0, then the operator 



Co 



is an exact symmetry of the equation 



F^{z) - 0. 



(2.10) 



(2.11) 



Remark 2.4. It is manifest from Eqs.( 2.8 ), (2.9 ) that if Xq is an exact symmetry of Eq.(2.11 ), then X = eXg 
is an approximate symmetry of Eq. (2.5). 



Definition 2.5. Eqs. (2.11) and (2.5) are termed an unperturbed equation and a, perturbed equation, respec- 
tively. Under the conditions of Theorem|2.3[ the operator Xq is called a stable symmetry of the unperturbed 



equation (2.11). The corresponding approximate symmetry generator X — Xq + eXi for the perturbed 
equation (2.10) is called a deformation of the infinitesimal symmetry Xq of Eq. (2.11) caused by the pertur- 
bation eFi{z). In particular, if the most general symmetry Lie algebra of Eq. (2.11) is stable, we say that 
the perturbed equation (2.5) inherits the symmetries of the unperturbed equation. 
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3 Approximate Symmetry of Gardner Equation 



Referring to Gardner equation ( |l.3[ ), if we put e = for small real parameter e, it becomes 

Wt- &{w + £w'^)Wx + Wxxx (3.1) 

Now, we can use Remar k 1 2 . 2 1 and Theorem|2.3|to provide an infinitesimal method for calculating approximate 



symmetries (2.6) for above differential equations with a small parameter. 



3.1 Exact Symmetries 

Let us consider the approximate group generators in the form 



d d d 

X = Xa + eXi = {aa+eai)— + (/3o + ePi)— + + £'?i)^, 

ox at aw 



where a^, j3i and rji for i = 0, 1 are unknown functions ol x,t and w. 
Solving the determining equation 

- 0, 



for the exact symmetries Xo of the unperturbed Gardner equation, means KdV equation we obtain 



where Ci,...,C4 are arbitrary constants. Hence, 



Vo = ^3 - 2C4W, 



00 O 

Xo - (Ci - 6C3t + Cix)— + (C2 + 3Cit)- + (C3 - ^C^w) — . 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



Therefore, unperturbed Gardner equation, means KdV equation, admits the four-dimensional Lie algebra 
with the basis 



^0^ 



dx' 



° ^ dx dw' 



d_ 
dx 



dt dw 



(3.6) 



3.2 Approximate Symmetries 



First we need to determine the auxiliary function H by virtue of Eqs.(2.8l, (2.9) and (2.5), i.e., by the 
equation 

'x^,''\F„{z)+eF,{z)) 



e 



Fo(z)+eFi{z)=o 



Substituting the expression (3.5) of the generator into above equation we obtain the auxiliary function 

H = l2wWx{C^w-C3). (3.7) 
Now, calculate the operators Xi by solving the inhomogeneous determining equation for deformations: 



H = 0. 



the above determining equation for this equation is written as 

XI '{Wt - 6wWa: + Wxxx) + l2wWx{CiW - C3) = 0. 

Solving this determining equation yields that C4 = 0, and hence, 

ai = C5 + Cgx - QCjt, Pi=Ce + iCgt, m ^ C7 - 2(^3 + Cs)^. 
Then, we obtain the following approximate symmetries of the Gardner equation : 



d d ^ d ,^ . d 

T^, V2 = — , V3 = 6t— + 2ew - 1 — , 
ox at ox aw 



(3.8) 



V4 = evi, V5 = ev2, vg = e(6i 



d__ 

dx dw ' 



ev3. 



V7 = eix 



d_ 
dx 



dt dw ' 



Because of C4 = 0, the scaling operator 

is not stable. Hence, the Gardner equation does not inherit the symmetries of the KdV equation. 

In the first-order of precision, We have the following Commutator table, shows that the operators (3.8) 
span an seven-dimensional approximate Lie algebra, and hence generate an seven-parameter approximate 
transformations group. 

Table 1: Approximate commutators of Approximate symmetry of Gardner equation 





Vl 


V2 


V3 


V4 


V5 


V6 


V7 


Vl 




















V4 


V2 








6vi 








6V4 


3v5 


V3 





— 6vi 








— 6v4 





-2V6 


V4 























V5 








6V4 














V6 





— 6V4 

















V7 


-V4 


-3v5 


2v6 















It is worth noting that the seven-dimensional approximate Lie algebra q is solvable and its finite sequence 
of ideals is as follows : 



C (V4) C (V4,V5) C (^V4,V5,V6^ C (^V4, V5, V6, V7 ) C ( Vi,V4,V5,V6,V7 
C (vi,V2,V4,V5,V6,V7) C 



4 Optimal System for Gardner Equation 

In general, to each s-parameter subgroup H of the full symmetry group G of a system of differential equations 
in p > s independent variables, there will correspond a family of group-invariant solutions. Since there are 
almost always an infinite number of such subgroups, it is not usually feasible to list all possible group- invariant 
solutions to the system. We need an efi'ective, systematic means of classifying these solutions, leading to an 
optimal system of group-invariant solutions from which every other such solution can be derived. 

Definition 4.1. Let G be a Lie group with Lie algebra g. An optimal system of sparameter subgroups is a 
list of conjugacy inequivalent sparameter subalgebras with the property that any other subgroup is conjugate 
to precisely one subgroup in the list. Similarly, a list of sparameter subalgebras forms an optimal system 
if every sparameter subalgebra of g is equivalent to a unique member of the list under some element of the 
adjoint representation: 1] = Ad{g{^)),g g G.[TT] 

Theorem 4.2. Let H and H be connected s-dimensional Lie subgroups of the Lie group G with corre- 
sponding Lie subalgebras [} and f) of the Lie algebra g of G. Then H = gHg~^ are conjugate subgroups if 
and only if f) = Ad{g{i))) are conjugate subalgebras. (Proposition 3.7 of 11 ) 



By theorem (4.2), the problem of finding an optimal system of subgroups is equivalent to that of finding 



an optimal system of subalgebras. To compute the adjoint representation, we use the Lie series 

2 

Ad{exp{fiVi))vj ^ Vj - /i[vj,Vj] + yK, [v^, v^]] - H = exp{ad{- ^iv i))v j , (4.1) 

where [vj, Vj] is the commutator for the Lie algebra, /i is a parameter, and i, j = 1, . . . , 7. In this manner, 
we construct the table with the (i,j)-th entry indicating Ad{exp{^'Vi)-v.j). 
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Table 2: Adjoint representation of Approximate symmetry of Gardner equation 



Ad 


Vl 


V2 




V3 


V4 


V5 


V6 


V7 


Vl 




Vl 


V2 




V3 


V4 


V5 


V6 


V7 — /iV4 


V2 




Vl 


V2 


V3 


- 6/iVi 


V4 


V5 


V6 - 6/iV4 


V7 - 3/iV5 


V3 




Vl 


V2 + 6/iVi 




V3 


V4 


V5 + 6/iV4 


V6 


V7 + 2/iV6 


V4 




Vl 


V2 




V3 


V4 


V5 


V6 


V7 


V5 




Vl 


V2 


V3 


— 6/iV4 


V4 


V5 


V6 


V7 


V6 




Vl 


V2 + 6/XV4 




V3 


V4 


V5 


V6 


V7 


V7 




Vl + MV4 


V2 + 3/XV5 


V3 


- 2/XV6 


V4 


V5 


V6 


V7 



Theorem 4.3. An optimal system of one-dimensional Approximate Lie algebras of the Gardner equation 
is provided by 

Vl, V4, V2+aV6, V3+Q!V2+/3V5 

vs + avi, Vfi + avi + /3v5, V7 + avi + /3v2 + 7V3. 

Proof. Consider the approximate symmetry algebra g of the Gardner equation, whose adjoint representation 
was determined in the table, our task is to simplify as many of the coefheients Ui as possible through judicious 
applications of adjoint maps to Vi. So that Vi is equivalent to V- under the adjoint representation. 
Given a non-zero vector 



Vl = aivi + 02 V2 + 03 V3 + 04 V4 + 05 V5 + aeve + a7V7. 

First suppose that 0,7 7^ 0. Scaling Vi if necessary, we can assume that a-j = 1. 
As for the 7th column of the table, we have : 

V'l = Ad{exp{{a4, - 3a6a5)vi))) o ylrf(ea;p(y V2)) o Ad{exp{-^v s))^i 

= (ai - 30206 - 2a3a5)vi + 02 V2 + 03 V3 + V7. 

The remaining approximate one-dimensional subalgebras are spanned by vectors of the above form with 
ar = 0. If as ^ 0, we have 

V2 = oivi -I- 02 V2 -I- V3 -I- 04 V4 -I- 05 V5 + aevg. 

Next we act on V2 to cancel the coefficients of vi, V4, vg as follows: 

= Ad{expQv2)) o Adiexpi^^^'^v,)) o Ad{exp{"^vr))V^ 

3a2a6 + 2a5 
= a2V2 + V3 H V5. 

if 03, 07 = and 02 ^ 0, the non-zero vector 

V3 = aiVi -I- V2 -I- 04 V4 -I- 05 V5 + CI6V6, 

is equivalent to: 

V3 = Ad{exp{-^V3)) o Ad{exp{ '^^"'^^^ ^^^^ ve)) o Ad{exp{^^vr))V3 = va + aeVg. 
if a2, as, ar = and ae ^ 0, we scale to make a^ = 1. Then 

V4 = aiVi -I- 04 V4 -I- 05 V5 -I- V6, 
is equivalent to V4 under the adjoint representation : 

V4 = ^rf(ea:p(^V2))V4 = aiVi + 05 V5 + vq. 
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if a2, as, «6, «7 = and 05 ^ 0, In the same way as before, the non-zero vector 

V5 = aivi + 04 V4 + V5, 

can be simplified : 

= ^d(exp(^V3))V5 = aiVi + V5. 
if a2,a3,a5,ae,ar = and ai 7^ 0, we act on 

Ve = vi + 04 V4, 

by Ad{exp{—a4V7)), to cancel the coefficient of V4, leading to 

Vg = ^rf(ea;p(-a4V7))V6 = vi. 

The last remaining case occurs when 01,02,03,05,06,07 = and 04 7^ 0, for which our earlier simplifica- 
tions were unnecessary. Because the only remaining vectors are the multiples of V4, on which the adjoint 
representation acts trivially. □ 



5 Approximately differential invariants for the Gardner equation 

In this section wc use two difi'erent methods to compute an approximately invariant solutions. 

In the beginning of this section we compute an approximately invariant solution based on the X = Vi +vq. 
The approximate invariants for X are determined by the equation : 



Equivalently : 



The first equation has two functionally independent solutions Jq = t and Jq=w. 

The simplest solutions of the second equation arc respectively, Ji = and Ji = x. Therefore we have 
two independent invariants w + ex = and ip{t) respect to X = Vi + Ve. 

Letting w + ex = ip{t) we obtain the w = ip{t) — ex for the approximately invariant solutions. 
Then, upon substituting in Gardner equation transform into the below equation : 

if' + 6e{ip — ex + e{ip — ex)^ = o{e). 

Therefore, in our approximation we have = ^p' = Q. Thus, invariant solution to Gardner equation Corre- 
sponding to X is : 

w = —ex. 

In this manner, we compute functionally approximate invariants respect to the generators of lie algebra and 
optimal system, as shown in Table 3 below. 

Where the unknown functions F, G, H are considered as follows: 

F{x, tM=xw( -l~2^+ + x^(-J^^ - , 
G{x,t,w) = , H{x,t,w) = — \-2{l -j)tw. 
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Table 3: Approximate Reduced equations corresponding to infinitesimal symmetries 



Operator 



Approximate Invariants 



Vl 
V2 
V3 
V4 
V5 
V6 
V7 

V2 + avg 

V5 + avi 
V3 + av2 + /3v5 
vg + avi + /3v5 

V7 + avi + /3v2 + 7V3 



t, w 

X, w 

t,x + 6tw — e{2xw + |j ) 
t, w 

X, w 

t,x + 6tw 

X — iat^e, w + aet 
at — ex, w if a 7^ 



x^w if a = 

— ax — e(3x, t + aw + £(— f — 2iw + Pw) 
( at — e/3x, aw + ex if a 7^ 
\ U"^ - (ix,t + fiw if a = 

^ 3tx 



a+e-yt^' 7(2; + 6tw) + aw + £F{x, t, w) 
at- (3x + 87*2 + eG{x, t, w),f3w + jt + eH{x, t,w) if /? 7^ 



ifa,/3,7 = 
if /3 = 0, + 7 V 



Unfortunately, for this method the first-order approximate generator does not necessarily yield a first- 
order approximate solution. The reason is that the dependent variables are not expanded in a perturbation 
series [T2]. 

Therefore we use another technique to find approximate invariant solutions for the Gardner equation. 
5.1 Approximate Galilean-invariant solution 

Now, we apply a different technique to find Approximate Galilean-invariant solutions for the Gardner equa- 
tion. We know that the general form of Galilean-invariant solutions to the unperturbed Gardner equation, 
means KdV equation, look as follows : 

W := 



6t 

a 



The function W is invariant under the operator 6tg^ — g^. 

we consider the approximate symmetry = V3 + fciV4 + fc2V5 -I- fcsvg + kiyj of the perturbed Gardner 
equation. Thus we will take the approximate symmetry 



d d { , d , d ^ , d d , , , d d d ,\ 

6t— - +s[ki— + k2Tr,+ h{Qt— -^) + ki{x— + 3t— - 2w— = Xo + eXi, 
ax ow V ox at ox aw ox at Ow / 

and use it for finding an approximately invariant solution by looking for the invariant perturbation of the : 



{53 



o(e). 



w ~W + ev{x, t). 
The invariant equation for above equation is : 

X^[w-W -ev{x,t)^ 

Note that Xq{w — W) vanishes identically. Therefore, Eq. ( |5.2[ ) becomes 

X^{w-W)~Xq{v)^Q, 
Therefore, we obtain the following differential equation for v{t, x) 



(5.1) 
(5.2) 



fcl + fc4C 



As you can see, the constant k^ is removed. It is easy to integrate this equation in the "natural" variables 

c — X 



W 



6t 



y = t. 
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Then it becomes : 

fci + kic + 6k2Z 

Vz H h 2z = 0. 

6y 

The integration yields 

(kiz + k4cz + 3k2Z^) o ^, , 

V = z +F{y). 

6y 

Returning to the variables t ,x, we have 

k2 ^'~f ~r-^Y+Fit). 



{ki + k4c){c ~ x) (c ~ x)'^ /r— T\2 



36^2 



Inserting this v in (5.1 ) and substituting in the perturbed Gardner equation we obtain t ^^^f^ + F{t) ~ 0, so 
Fit) = f , where C is an arbitrary constant. 

Thus, the approximate symmetry provides the following the approximately Galilean-invariant solu- 
tions : 



uj^W + e{ ik,+k,c){c-x) _ jc-xf fc-x^\C^ 



(^) +t)- 



Galilean-invariant solution, means W , and Approximate Galilean- invariant solution, means w, are diplayed 
below for — 3 < a; ^ 3, 0.1 ^ t ^ 3, c = C = /ci = fc2 = A:4 = 1 and e — 0.1, respectively. 





(a) Galilean-invariant solution (b) Approximate Galilean-invariant solution 

Figure 1: Comparison of Galilean- invariant and Approximate Galilean-invariant solution 
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